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ABSTRACT 

Wc examine the effects of gas expulsion on initially sub-structured and out-of- 
equilibrium star clusters. We perform iV-body simulations of the evolution of star 
clusters in a static background potential before removing that potential to model gas 
expulsion. We find that the initial star formation efficiency is not a good measure of 
the survivabihty of star clusters. This is because the stellar distribution can change 
significantly, causing a large change in the relative importance of the stellar and gas 
potentials. We find that the initial stellar distribution and velocity dispersion are far 
more important parameters than the initial star formation efficiency, and that clusters 
with very low star formation efficiencies can survive gas expulsion. We suggest that it 
is variations in cluster initial conditions rather than in their star formation efficiencies 
that cause some clusters to be destroyed while a few survive. 
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1 INTRODUCTION 

Many, perhaps the vast majority, of stars form in dense, 
clust ered environments (|Lada fc Ladal [20031 : iBressert et al.l 
l20ld ). However, after only 10 - 20 Myr the vast majority of 
stars are dispersed into t he low density environment of the 
field (jLada fc Ladall2003h . The destruction of most of these 
early dense clusters has been termed 'infant mortality'. 

Perhaps the best candidate for the mechanism behind 
infant mortality is primordial gas expulsion. Giant molec- 
ular clouds turn only a few per cent to a few t ens of per 
cent of their gas into stars (jLada fc Ladal l2003l ) meaning 
that the potentials of very young embedded clusters 
are dominated by gas. Feedback from the most massive 
stars will remove this gas on a timescale of a few Myr 
significantly altering the potential of the cluster and 
potentially leading to its destruction. This process has been 



ceptions are iLada et al.lll984l: lElmegreen fc ClemendllQSSl : 



IVerschueren fc Davidll989l : lGoodwinll2009l ) and structurally 
simple (usually a Plummer sphere or similar). But re- 
cent observational and theoretical results strongly suggest 
that stars do not form in d ynamical equilibrium, nor in a 



smooth distribution (see e.g. Elmegreen fc ElmegreerJ 2001 
Allen et all l2007l : iGutermuth et all 120091 "! Bressert et al.l 



2010l and references in all of these papers) . This should not 
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However, previous work has tended to concentrate on 
clusters that are both in virial equilibrium (notable ex- 
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be surprising as in the gravoturbulent model of star forma- 
tion, stars will form in dense gas in filaments and clumps i n 
a turbulent environment (see e.g. [Elmegreen fc Scalci|[2004l ). 
This has lead to a model of star cluster formation in which 
a dynamically cool and clumpy initial state collapses and 
violently relaxes into a d ense star cluster (e.g. lAllison et al.l 
|2009| . lAllison et ai]|2010l ') . 

Strongly non-equilibrium initial conditions may play a 
crucial role in the effects of gas expulsion on clusters be- 
cause the key factor in determining the effects of gas expul- 
sion is not only the depth of the gas potential, but also the 
dynamical state of the stars at the onset of gas expulsion 
IVerschueren fc David|[l989l : lGoodwin|[2009l l . As an extreme 
example, if 99 per cent of the mass of a cluster is expelled, 
the cluster will still survive if the stars have almost no ve- 
locity dispersion. 

We might therefore expect the effects of gas expulsion 
on a cluster to depend on three factors. Firstly, a clumpy 
and non-symmetrical initial distribution of the stars. Sec- 
ondly, the star formation efficiency (SEE), and so the rel- 
ative importance of the stellar and gas potentials. Thridly, 
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the (non-equilibrium) initial conditions of the cluster and 
the evolution of these initial conditions up until the point 
of gas expulsion. If clusters are unable to relax fully before 
gas expulsion then we might expect them to respond very 
differently to the loss of the residual gas. 

In this paper we examine the A^-body evolution of 
highly non-equilibrium star cluster in fixed background po- 
tentials to model the background gas. The background po- 
tential is then removed to simulate the effects of gas ex- 
pulsion. We describe our initial conditions in Section 2, our 
results in Section 3, and then discuss the potential conse- 
quences in Section 4 before drawing our conclusions in Sec- 
tion 5. 



2 INITIAL CONDITIONS 

We p erform our A ^-body simulations using the Nbody2 
code (|Aarsethll200ll V Nbody2 is a fast and accurate direct- 
integration code, optimised for the number of star particles 
used in this study (A*' = 1000). In our simulations there are 
two separate mass distributions which we wish to model: the 
stars and the background gas potential. 



2.1 Stellar distribution 

In all cases we model the stellar distribution as N = 1000 
particles with equal masses of O.5A/0 resulting in a total 
stellar mass of 500^/0. Every particle has a gravitational 
softening length of 100 AU. We choose softened and equal- 
mass particles in orde r to avoid strong two -bod y interactions 
and m ass segregation. lAUison et al.l (|2009l ) and lAUison et all 
(|2010D showed that both of these effects can be extremely 
important in the violent collapse of cool, clumpy regions, 
however we wish to avoid complicating our simulations with 
these effects as we are interested in the effects of gas expul- 
sion. 

We distribute the stars within a radius of 1.5 pc with 
two different 'clumpy' morphologies: 'fractal', and 'clumpy 
plummer'. Representative snapshots of the initial conditions 
of a fractal and clumpy plummer distribution can be seen in 
the left and right-hand panels, respectively, of Fig. (2] 

Fractal cluster s are produced using the box fr actal tech- 
ni nue described bvlGoodwin fc Whitwortl] (|2004| ): also used 
bv Ullison et all (|2009l ') and lAUison et all ((2OI0I ). In this pa- 
per we use a fractal dimension D = 1.6 which corresponds 
to a highly clumpy initial distribution (see Fig. [2]). 

In clumpy plummers, each sub-clump is modelled as 
an individual Plummer sphere. The total stellar mass of 500 
M0 is distributed into 16 sub-clumps of approximately equal 
mass. Each sub-clump contains 31.0-31.5 Mq formed from 
62-63 particles. The position of individual sub-clumps within 
the model star forming region follows the gas potential (see 
below) . 



2.2 Gas distribution 

The gas within the star-forming region is modelled as a static 
Plummer potential with a mass Mg and scale-radius rg. 

The scale radius is set to be rg = 1 pc, or rg = 1.5 pc. 
As the stars have a maximum radius of 1.5 pc this means 



that for rg = 1.5 pc, the gas distribution seen by the stars 
is roughly uniform. 

The relative masses of the gas potential within 1.5 pc 
and the stars set the true SEE (e), that is the fraction of the 
initial cloud mass converted into stars. The mass of the gas 
potential is varied to obtain true SFEs of 20, 30 and 40 per 
cent (Mg =2000, 1167, 750 Mq). 

Note that a true SEE of 20 or 30 per cent will fail to 
produce a bound star cluster after instantaneous gas ex- 
pulsion if the gas and stars are initially virialised (see e.g. 
iBaumgardt fc Kroupall2007l ). 

We emphasise that the gas potential does not follow 
the initial stellar distribution, nor does it react to changes 
in the stellar distribution as it evolves (it is not live). These 
are obviously extreme simplifications, but as we will discuss 
later we feel that we capture the essence of the basic physics 
using such a simple model. 



2.3 The dynamical state of the stars 

The virial ratio Q — T/\il\ is the ratio of the kinetic , T, to 
potential energy, Q,, of the system. We set the initial velocity 
dispersion of the stars relative to the total potential (gas & 
stars) to set initial stellar virial ratios between zero (cold), 
0.5 (virial equilibrium), and 0.9 (supervirial, but bound). 

For initial virial ratios Q < 0.5, the system will tend to 
collapse, and for Q > 0.5 the system will tend to expand. 
However, even for Q — 0.5, although the system is in virial 
equilibrium, the clumpy initial conditions mean that it is 
not in dynamical equilibrium. 



2.4 Gas expulsion 

Although young star clusters form embedded within the 
the molecular gas from which they formed, few star clus- 
ters over ~ 5 Myr old r emain associated with their gas 
(iProszkow fc Adamd [20091 ). This is likely as a result of a 
number of mechanisms including radiative feedback from 
massive stars, stellar winds from young stars, and eventually 
the onset of the first supernova. The time at which gas re- 
moval begins to occur, and the duration of the gas removal 
process is uncertain, and dependent on the particular gas 
removal mechanism in operation. 

To simulate gas expulsion we instantaneously remove 
the external gas potential after 3 Myr. This is approximately 
mid-way between the time at which gas removal from stellar 
winds and supernova feedback might be expected to occur 
for star clusters containing ~ 1000 stars. Instantaneous gas 
removal is the most extreme form of gas removal as the stars 
have no chance to readjust to the change in the potential 
re.g. lGood^ll997al : [Baumgardt fc Kroupall2007l ). 



2.5 Ensembles 

It should be noted that both the clumpy plummer and frac- 
tal clusters can vary considerably in appearance depending 
on the random realisation u sed and the subseq uent evolu- 
tion is highly stochastic (see I Allison et aLll2010l ). We there- 
fore conduct a minimum of 5 random realisations of each 
parameter set. 
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Figure 1. An x-y projection of the initial (left panel) and final (right panel) stellar distribution of an initially clumpy plummer morphology 
stellar distribution, evolved for 3 Myr in an rg = 1 pc, Mg = 2000Mq gas potential. The initial virial ratio of the stellar distribution is 
0.35. Time of snapshot and LSF at that instant is indicated in top-right hand corner of each panel. The dotted (red) circle marks the 
stellar half-mass radius, directly used to calculate the LSF. In this case, the LSF almost doubles from 0.23 (left panel) to 0.40 (right 
panel) over the duration of the embedded phase. 



Total mass, Mtot 


1250-2500 Mq* 


Outer radius 


1.5 pc 


Star formation efficiency, e 


20-40%* 


Total stellar mass, Mstar 


500 M0 


Particle number, N 


1000 


Particle mass, rripart 


0.5 Mq 


Stellar morphology 


fractal or Plummer * 


Initial virial ratio, Q; 


0.0-1.0* 


Total gas mass, Mg 


750-2000 M0* 


Gas Plummer scale radius, r-g 


1.0 pc or 1.5 pc* 


crossing-time, Tct 


1.4 Myr/2.0 Myr 




(rg=1.0 pc/1.5 pc) 



Table 1. The initial conditions parameter sets. Parameters are 
loosely divided into quantities governing; the total region (up- 
per), the stellar component (middle), the gas component (lower). 
Parameters that arc highlighted with a '*' are parameters we ad- 
just between simulations. 



2.6 Summary 

We set-up clumpy SOOM© star clusters with A'' = 1000 equal- 
mEiss particles within a static background gas potential. The 
dynamical state of the stars varies from very cold to al- 
most unbound. The stars dynamically evolve within the gas 
potential before its instantaneous removal after 3 Myr. We 
measure the final properties of each star cluster an addi- 
tional 3 Myr after gas expulsion. A summary of the key 
parameters is provided in Table [T] 



3 RESULTS 

We will first discuss the evolution and relaxation of the stel- 
lar distribution in the gas potential before addressing the 
response of the stars to gas expulsion. 

3.1 The embedded phase 

In all simulations, we evolve our stellar initial conditions 
for 3 Myr in the gas potential. In Figure |3] we present an 
example of the evolution of an initially fractal stellar distri- 
bution with an initial virial ratio of 0.33 in a gas potential 
of scale-radius 1 pc. The initial true SFE of this system is 
e = 20%. 

Clumpy stellar sub-structure interacts resulting in 
clump merging, and stellar scattering from the clumps 
(somewhat suppressed in these simulations due to the soft- 
ening). As a result, substructure is significantly erased on 
~l-2 crossing-times of the region. For example, the crossing- 
time of the region presented in Fig. |3]is 1.4 Myr. By 3 Myr 
(bottom-right panel) over 2 crossing-times have elapsed, 
and the remaining stellar distribution has lost most of its 
su b-structure. Similar behaviour i s seen in the s imula tions 
of 'Goodwin & Whi tworthI (|2004l ). lAUison et ahl (|2009l ) and 
lAUison ct al. (2010J. However, they do not include a model 
for the gas component of the region in their simulations. By 
including a gas potential in our simulations, the crossing- 
time of the region is shortened. However the picture of evo- 
lution from clumpy substructure to an almost smooth distri- 
bution of stars, on timescales of ~l-2 crossing-times of the 
region, is unsurprisingly unchanged by the presence of the 
gas. 

Importantly, the relaxation and smoothing of the stellar 
distribution has resulting in it shrinking significantly (see 
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t = 0.0 Myr 




Figure 2. Representative examples of an initial distribution of 
morphology (right panel). 

also lAllison et al] l2009l : I Allison et al]|2010l '). It is clear in 
Fig. [3] that the central regions of the cluster now contain a 
greater number of stars than initially. This means that (be- 
cause the gas potential is static) that the stars are a more sig- 
nificant contributor to the cent ral potential than initially - 
the effective SF E has increased jVerschueren fc David|[l989l : 
lGoodwin|[2009l ) and, in this case, the cluster should be more 
able to survive gas expulsion. 

In order to quantify the effects of collapse (or expansion) 
and the change of the effective SFE we define the Local 
Stellar Fraction (LSF) of the star forming region: 

LSF = A/^(r < nO 

A'L{r < rh) + Mg^^r < n,) ' 

where rh is the radius from the centre of the region contain- 
ing half the total mass of stars. M* and Mgas is the mass of 
stars and gas, respectively, measured within rh. Therefore 
the LSF is a measure of the current effective SFE. 

In any non-equilibrium system the half-mass radius of 
the stars will change with time. Figure [T] shows a projection 
of the stellar distributions of an initially clumpy plummer 
model within an rg = 1 pc, Mg = 2OOOAf0 gas potential. 
Initially the stars have an initial virial ratio of 0.35, and a 
half-mass radius r-h — 0.78 pc. LSF — 0.23, close to the 
initial true SFE of e = 0.2 (left panel). However, at 3 Myr 
the stellar half-mass radius has fallen to rh ~ 0.55 pc and 
the LSF is 0.73 (right panel) - almost double the initial 
value. 

In Figure [4] we plot the LSF measured moments be- 
fore gas expulsion, versus the initial virial ratio of the stars 
Qi. Unsurprisingly, stellar distributions with low initial virial 
ratio form clusters with high LSFs. What is slightly surpris- 
ing is that high-Qi (warm) clusters do not have significantly 
lower LSFs than their initial true SFE. This is because the 
gas potential dominates and even almost unbound clusters 
only expand by a factor of less than two in their half-mass 
radii (some stars are lost, but the half- mass radii do not 
change by very significant factors). However, there is con- 
siderable scatter in the figure. The values with errorbars 



t = 0.0 Myr 

1 - 




with a fractal morphology (left panel), and a clumpy plummer 



connected by a solid line in Figure |4] represent the mean 
values with standard deviation of this scatter in LSF (for 
e = 0.2 simulations) measured in 0.1 wide bins of Qi. We 
show a representative histogram for the (5i=0.3-0.4 bin in 
Figure [5] On top of our 5 random realisations of each point 
of our parameter space, we conduct an additional 60 Plum- 
mer and 70 fractal morphology simulations so as to ensure 
a minimum of 20 simulations in each Qi bin to avoid low 
number statistics. The scatter is higher for stellar distribu- 
tions with low initial virial ratios. This scatter reflects the 
stochastic nat ure of cluster evolu tion from clumpy initial 
conditions (see lAllison et alll2010l ). 

We find no evidence for any dependency in the position 
or scatter of the trend on initial stellar morphology (clumpy 
plummer or fractal) or gas concentration (uniform or con- 
centrated). There is a trend in that clusters with higher true 
SFEs have higher LSFs reflecting their initially higher SFEs 
(open boxes and triangles in Figure |4]). 

As well as the LSF, the virial ratio of the stars at the 
onset of gas expulsion is also important. In particular, the 
velocity dispersion of the stars is crucial, with sub- virial 
clust ers much more able to survive gas expulsion (jGoodwinl 
I2OO9I ). Therefore we measure the stellar virial ratio just be- 
fore gas expulsion Qpc- 

Figure[5]shows the initial and pre-expulsion virial ratios 
for all our clusters. As would be expected, our clusters have 
attempted to reach equilibrium (by erasing substructure and 
collapsing/expanding). In all cases the clusters have reached 
an approximate virial equilibrium of Qpc ~ 0.5. However, it 
is crucial to note that these are approximate equilibriums, 
and none of these clusters are fully virialised. We shall return 
to this point later. 



3.2 Post-gas-expulsion evolution 

After 3 Myr of evolution within the embedded phase (i.e. 
within the gas background potential), we assume the gas is 
instantaneously expelled from the star forming region. This 
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Figure 3. The evolution of an initially fractal {D = 1.6) stellar distribution in an Tg = 1 pc, Mg = 2000Mq gas potential. The initial 
virial ratio of the stellar distribution is 0.33. The time of each snapshot is indicated in top-right hand corner of each panel. 
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0.2 0.3 0.4 0.5 0.6 0.7 O.S 

Initial virial ratio Q. 

Figure 4. Plot of the initial virial ratio Q; against the Local Stel- 
lar Fraction [LSF] at 3 Myr. Symbols reflect the initial true SFE 
of the region of 20, 30 or 40 per cent. Initially hot dynamical tem- 
peratures result in low LSF, and conversely very cold dynamical 
temperatures result in high LSF. At a fixed value of initial virial 
ratio, there is a slight trend for raised LSF for increasing true 
SFEs e. Error bars show the mean and standard deviation of the 
scatter observed in the simulations with an initial true SFE of 20 
per cent. 




Figure 5. A representative histogram of the scatter in LSF within 
Qi bin 0.1-0.2. Solid bars represent the numbers counted within 
each LSF bin. The cross point is the average LSF for all clus- 
ters within the bin. The error bars of the crossed point arc the 
standard deviation of the LSF scatter. The dashed line shows a 
gaussian fit with equal standard deviation. 



is simply modelled by an instantaneous removal of the gas 
potential. We further evolve the gas-free star cluster an ad- 
ditional 3 Myr after gas expulsion. At this time, we measure 
the bound mass of the star cluster. Any star whose kinetic 
energy is less than it's potential energy is considered bound. 
We normalise the bound mass by the total stellar mass in 
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Figure 6. Comparison of the initial virial ratio Qi, chosen at the 
beginning of the simulation, versus the pre-gas-expulsion virial 
ratio Qpe, measured just prior to gas expulsion. Pre-expulsion 
virial ratios are all close to virialised (Qpe = 0.5) but scatter 
around this value with a standard deviation ~ 0.05. 



the simulation - this quantity is referred to as the bound 
stellar fraction (./bound)- 

In Figure [7] we plot the bound stellar fraction /bound, 
as a function of t he initial true SFE e. We also plot the 
results of .Baumgar dt fc Kroupal (|2007l ') (herein BK07) who 
use initial smooth and virialised equilibrium clusters as their 
initial conditions. 

It is quite clear from this figure that our clusters with an 
initial true SFE of 20 or 30 per cent can survive gas expulsion 
with very significant fractions of their initial mass remaining. 
Also there is a very significant scatter in our results rather 
than the very tight relationship found by BK07. 

We confirm that our choice of 3 Myr of gas free evo- 
lution has not significantly infiuenced our measured bound 
fraction by evolving a small sub-sample of clusters for 15 
Myr after instantaneous gas removal. We continue the gas 
free evolution phase of three clusters that result in a high 
(/bound = 0.72)), medium (/bound = 0.45)) and low mass 
cluster (/bound = 0.30)) when measured after 3 Myr. After 
15 Myr, the measured bound fraction differs by < 5% from 
those measured after 3 Myr. 

The differences in our results from those of BK07 are in 
fact due to the differences in the initial conditions used for 
the simulations. As we have seen in the previous section, the 
effective SFE as measured by the LSF can increase signif- 
icantly before gas expulsion. Rather than expecting a rela- 
tionship with the initial SFE we should expect a relationship 
with the LSF which we plot in Fig. [8] Here we find that the 
final bound fraction of stars does correlate with the LSF in 
a similar way to the results of BK07. Note that there is no 
systematic differences between the results of simulations of 
fractals or clumpy plummers, nor with the concentration of 
the gas potential. 

However, Figure [8] shows there is still a considerable 
scatter in the fhound-LSF trend. Indeed, some clusters with 
an LSF of < 0.3 can retain a significant bound fraction of 
stars after gas expulsion (up to 50 per cent). This is because 
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Star formation efficiency e 

Figure 7. Plot of the region's SFE e versus the final cluster's 
bound mass fraction /bound- Crossed symbols are the results of 
all simulations in this paper, where as the solid line is the results 
of the BK07 simulations. The points are highly scattered with 
no clear trend, and do not match the results seen in the BK07 
simulations. 
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Figure 9. Plot of the influence of the dynamical state of the 
cluster (as measured by Qpe) at the time of gas removal. All 
symbols in Figure \8\ (the standard model simulations) are now 
shown as crosses. Circle and square symbols represent simulations 
where Qpe > 0.55 or Qpe < 0.45 respectively as indicated in the 
key. 
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Figure 8. Plot of the influence of LSF (measured at the onset of 
gas expulsion) on the flnal bound mass fraction (/bound) of clus- 
ters. All open symbols represent standard model simulations with 
varying initial stellar morphology or gas potential concentration. 
Key labels are fractal (f) or clumpy plummer (p) morphology, 
and uniform (u) or concentrated (c) gas concentration. There is 
a good trend for increasing bound mass fraction with increasing 
LSF, although clearly there is scatter in the trend. The trend is 
independent of initial stellar morphology or gas concentration 



the LSF is not the only parameter influencing stellar mass 
loss as a result of gas expulsion. The dynamical state of the 
cluster, as meas ured by Qpc, is also expected to be important 
(|Goodwinll2009l ). Recall from Figure [6]- by the onset of gas 
expulsion all clusters had become approximately virialised. 
However, even minor deviations from virial equilibrium at 
the onset of gas expulsion can significantly effect the cluster 



evolution post-gas-expulsion. In Figure |9l we again plot the 
fhound-LSF trend, but this time highlighting simulations 
whose pre-gas-expulsion virial ratio Qpc is greater than 0.55 
(super- virial) or less than 0.45 (sub- virial) . Clusters that are 
slightly super- virial at the onset of gas expulsion clearly lie 
to the lower bounds of /bound, whilst those which are sub- 
virial are at the higher-end. 

A gravitational system will reach approximate virial 
equilibrium in 1 - 2 crossing times, but it will tend to oscil- 
late around exa ct virial equilibrium for some time (see also 
lGoodwin|[T997al ). In Figure [TOl we plot the evolution of the 
virial ratio with time during the embedded phase. We in- 
clude representative curves for simulations with initial virial 
ratios of Qi = 0.0, Qi = 0.33, Qi = 0.47 and Qi = 0.6. 
All curves are results from simulations with a concentrated 
gas profile with a crossing-time of 1.4 Myr. Whilst all clus- 
ters reach approximate virial equilibrium quite quickly they 
oscillate around Q = 0.5 for some time. And, as shown in 
fig. [9l the virial ratio they have at the moment of gas ex- 
pulsion can make a significant difference to their ability to 
form a bound cluster afterwards. 



4 DISCUSSION 

We have shown that initially clumpy and out-of-equilibrium 
clusters will rapidly relax and erase most of their substruc- 
ture. Depending on their initial virial ratios and substruc- 
ture this may mean that they may either collapse or expand. 
This will alter the importance of the gas potential on the 
stars meaning that initially cool clusters are more able to 
survive gas expulsion, whilst initially warm clusters are less 
likely to. This can remove much of the dependence of sur- 
vival after gas expulsion from the importance of the initial 
true SFE and place it instead on the i mportance of the initial 
dynamical state of the stars (see also IVerschueren fc DavidI 
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Time / Myr 



Figure 10. Time evolution of tiie virial ratio for Qi = 0.0 (solid, 
black lines), Qi = 0.33 (dot-dash, blue lines), Qi = 0.47 (dotted, 
green lines), and Q; = 0.6 (dashed, red lines). Colour version 
available online 



1 19891 : lGood^l2009l ). We have also shown that there is a 
very significant scatter due to both the intrin sic differences 
betw een (statistically the same) clusters (see lAUison et al.l 
I2OIOI ). and on the exact virial ratio at the onset of gas ex- 
pulsion. 

There is significant observational and theoretical evi- 
dence that the initial distrib utions of stars are not smooth , 
nor are th e y viri a lised (see Elmegreen fc Elmegreen 200 ll: 
Bate et al.l I2OO3I: iBonnell et al.l l2003l: I Bert out fc Genoval 



20061: lAUenet al.1 l2007l : iKraus fc Hilleribrandl I2OO8I 



Gutermuth et al. 
2OIOI and references in 



I2OO9I : IClarkd l2oTo[ iBressert et all 



natu ral consequence of 



e.g. lElmegreen fc ScaloT l2004l: 



all of these papers) which is a 
ravoturbulent star formation (see 



McKee fc Ostrikeij l2007l : 



2OIOI ). Indeed, evidence 



iBergin fc Tafallal |2007|~ IClarkd 

points towar d s sub -virial initial conditions for stars (see 
ison et all (j2010l ) and references therein). We therefore 
argue that our initial conditions are far more realistic 
than those o f a smooth, relax e d cluster as generall y used 



before (e.g. [ Goodwin' 'l997al: iGoodwin fc BastianI l2006l : 



iBaumgardt fc Kroupai2007l). 

iProszkow fc Adamd (|2009l ) conduct a large parameter 



study investigating the key parameters controlling the final 
bound fraction of clusters that have undergone gas expulsion 
using smooth and spherical initial stellar distributions. They 
find a clear trend with SFE although they, too, use sub-virial 
initial stellar distributions. We note that we also see a trend 
for increasing bound fraction with increasing SFE (Figure 
O, but that it is highly scattered. The source of this scatter 
is the use of clumpy initial stellar distribution. Therefore the 
strength of the star formation efficiency as a predictor for 
the survival of a cluster to mass loss is severely weakened 
with the use of far more realistic initially clumpy stellar 
distributions. 

It has often previously been assumed that for a cluster 
to survive it must have had a high SFE and so the small 
numbers of clusters we s ee must be a high-SFE tail to the 
SFE distribution (see e.g. IParmentier et al.ll200^ ). However, 



we have shown that some low-SFE clusters can survive if 
they are 'lucky' enough to have the right initial conditions. 
Indeed, the low sur vival rates found fo r young clusters of 
only ~ 10 per cent (|Lada fc Lada|[2003h may be better ex- 
plained as these being the few clusters with the 'right' initial 
conditions than being an extremely high-SFE (> 40 or 50 
per cent) tail of star formation. 

We note that our si mulations are h ighly idealised. Stars 
are equal mass, when lAUison et al.l (|2009l ) have demon- 
strated that rapid mass segregation can occur with a more 
realistic initial mass function. 

We do not include any primordial binaries, although 
their presence can clearl y influence the dynam ics of a cluster 
( Goodman fc Hull Il989l ). The simulations of iKroupa et al.l 
( 20011 ) show that changes in binary fraction, and scattering 
between stars in the stellar outflow (that can result in bi- 
nary hard ening) can increa s e the flnal bound fraction of a 
cluster. In lMoeckel fc Bat3 (|2010l ) the binary fraction does 
not change signiflcantly during the cluster formation process 
as a result of formation in an initially high stellar density 
environment - within sub clumps. However their initial con- 
ditions are limited to statistics of one, and we have further 
demonstrated that clumpy initial conditions can result in 
highly stochastic behaviour. 

Furthermore we assume in s tanta neous gas removal, al- 
though [Baum gar dt fc Kroupal l|2007h demonstrate that a 
slower rate of gas removal can result in a higher cluster sur- 
vival rate. 

The length of the embedded phase is flxed at 3 Myr in 
our simulations although this could vary depending on the 
nature of the gas removal mechanism. Despite these simplifi- 
cations, we argue that the idealised nature of the simulations 
has enabled us to more clearly test the implications of an 
initially sub-virial, and clumpy stellar distribution. We defer 
a less idealised study to a later paper. 

Our simulations do have an obvious problem, however, 
in that we use a smooth, static background potential for the 
gas. This has two main problems. Firstly, the initial gas and 
stellar distribution do not match, despite the fact that our 
stars are assumed to have formed from this background gas. 
Secondly, the gas is not able to respond to the motion of 
the stars. (A third, but less important problem, is that we 
assume that all of the stars form instantaneously.) 

The importance of both problems comes down to how 
well the motions of the gas and the stars are coupled. If 
conditions are such that both stars and (at least a significant 
fraction of the) gas move together in the potential then we 
would expect both to collapse or expand together and the 
LSF and true star formation efficiency to remain roughly 
constant. However, if the bulk of the gas does not notice the 
stars because it is not involved in their formation and the 
relative gravitational influence of the stars is small, then 
our approximations should be roughly correct. We would 
argue that at low true star formation efficiencies that the 
bulk of the gas would be uncoupled from the stars. We are 
working on more detailed simulations with a live background 
potential which we will present in future papers. 
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5 SUMMARY & CONCLUSIONS 

We perform A''-body simulations of sub-structured, non- 
equilibrium 5OOM0 clusters of A'^ = 1000 equal-mass stars 
in a static background potential. The mass of gas is varied 
to simulate star formation efficiencies (SFEs) of 20 to 40 per 
cent. After 3 Myr of dynamical evolution, the potential is 
instantaneously removed to model the effect of gas expulsion 
from the cluster. 

Previous work with initially smooth and equilibrium 
clusters has shown that there is a critical SFE for the sur- 
vival of (at least par t of) t h e cluster of ~ 30 per ce n t (e.g . 
iGoodwin fc BastianI (|2006l ): iBaumgardt fc Kroupal (|2007l ) 
and references therein). However, it has also been known 
that it is the conditions at the onset of gas expulsion that are 
important in influenci ng the evolution of the star cluster fol- 
lowin g gas expulsion jVerschueren fc David|[l989l . iGoodwinI 
l2009l ). 

Our key results may be summarised as follows. 

(i) The true SFE is not a good measure of the survival 
or otherwise of initially out-of-equilibrium clusters. Even 
though clusters rapidly come into approximate virial equi- 
librium, their structure can change significantly and, unless 
the gas follows the motions of the stars, the ratio of gas 
mass-to-stellar mass quantified by the local stellar fraction 
(LSF) can change hugely. 

(ii) It is the LSF, measured at the onset of gas expul- 
sion, that is the key parameter controlling the survival of an 
embedded cluster to gas expulsion (see Figure [8|. A cluster 
with a high LSF at the onset of gas expulsion - regardless 
of the initial SFE - is able to produce a cluster with a signif- 
icant fraction of the initial cluster mass after gas expulsion. 

(iii) The dynamical state of the cluster, as measured by 
the virial ratio at the onset of gas expulsion, is also im- 
portant. If gas expulsion occurs when the cluster is mildly 
sub- virial (i.e. collapsing), the cluster will suffer lower stellar 
mass loss as a result of gas expulsion. Conversely, a cluster 
that is mildly super-virial (i.e. expanding) at the onset of 
gas expulsion suffers greater stellar mass loss during gas ex- 
pulsion. 

The initial SFE of a cluster is almost certainly not a 
good measure of the ability of a cluster to survive gas ex- 
pulsion. The initial spacial and kinematic distributions of 
the stars is crucial, as is how the stars and gas both evolve 
once the stars have formed. It is quite possible for low-SFE 
clusters (< 20 per cent) to produce bound clusters after gas 
expulsion given the right initial conditions. 
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